THE EHRHART FUNCTION FOR SYMBOLS 
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Abstract. We derive an Ehrhart function for symbols from the Euler-MacLaurin 
formula with remainder. 
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1. Introduction. 

Let A C R™ be a convex polytope whose vertices are in Z" and such that the 
origin is in the interior of A. Consider the expanded polytope N ■ A. 

Ehrhart's theorem |Ehr| asserts that for N a positive integer, the number of 
lattice points in the expanded polytope, i.e. 

#((iV-A)nZ"), N G Z + 

is a polynomial in N. More generally, suppose that / is a polynomial, and let 



(1) 



p(N, f) 



Then Ehrhart's theorem asserts that p(N, /) is a polynomial in N. 

In the case that A is a simple polytope (meaning that n edges emanate from 
each vertex) Ehrhart's theorem is a consequence of the Euler-MacLaurin formula, 
|Kb1 iKPl RlSn IHS21 IHul IPTviriTR] and one can be more explicit about the nature 
of the polynomial p(N, f). 

Let us explain how this works in the more restrictive case where A is not only 
simple but is regular, meaning that the local cone at each vertex can be transformed 
by an integral unimodular affine transformation into a neighborhood of the origin 
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in the standard orthant R™ - 1 In this case we can apply the formula of Khovanskii- 
Pukhlikov |KPj . which reads as follows: The polytope A can be described by a set 
of inequalities 

x ■ Ui + cij > 0, i = 1, . . . , m 

where m is the number of facets of A, where the Ui are primitive lattice vectors, and 
where the di are positive integers. Then for any positive number t, the expanded 
polytope t ■ A is described by the inequalities 

x ■ Ui + tai > 0, i = 1, . . . , m. 

Let A tj /j, h = (h\, . . . , h m ) be the polytope defined by 

(2) x ■ Ui + tai + hi > 0, i = 1, . . . , m. 
Then the function 

(3) P(t,h,f) := f f(x)dx 

JAt.h 

is a polynomial in t and h. 

The formula of Khovanskii-Pukhlikov (applied to N ■ A) expresses p(N, /) in 
terms of a differential operator applied to p(t, h, /): Explicitly, consider the infinite 
order constant coefficient differential operator 

< 4 > ™d(£)=v>(^ 7 

where 6 7 x 7 is the Taylor series expansion at the origin of the Todd function 

m 

i=i 

The Khovanskii-Pukhlikov formula says that 

p(NJ)= Todd (Jj^p(N,h,f) ^ . 

Note that since p is a polynomial in h the right hand side really involves only a 
finite order differential operator. 

For purposes below it will be convenient to write the Khovanskii-Pukhlikov for- 
mula in the form 

(5) p(N, /) - p(N, 0, /) = (^Todd (J^ - Id) p(N, h,f) . 

For simple polytopes there is a more general formula due to |CS2llUullBV| . Our 
goal in this paper is to prove an analogue of |5J and its generalizations when the 
polynomial / is replaced by a "symbol" , a term whose definition we recall from the 
theory of partial differential equations: A smooth function / S C°°(R™) is called a 
symbol of order r if for every n-tuple of non-negative integers a := (ai, . . . , a„), 
there exists a constant C a such that 

\d?...d$rf(x)\<c a o. + \x\y-w 



^Of course, this use of the word "regular" has nothing to do with the term denoting Pla- 
tonic solids. There are other names in the literature for the property we are describing, such as 
"smooth", "Delzant" , "torsionfrec" , "unimodular" etc.. It is unfortunate that the nomenclature 
for polytopes with this property has not yet been standardized. 
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where \a\ = ai. In particular, a polynomial of degree r is a symbol of order r. 
Note that if / is a symbol of order r on R™ then its derivatives of order a are in L 1 
if r < \a\ — n. 

For simplicity, we will restrict ourselves in this paper to polyhomogeneous 
symbols, meaning functions / G C°°(R™) which admit asymptotic expansions of 
the form: 

r 

(6) /(*) 

— oo 

for 1 1 a; 1 1 >> where the fi are homogeneous symbols of degree £ meaning that each 
fl is a symbol of order I with the property that for ||x|| sufficiently large and t > 1 
we have ftitx) = t l fi{x). The sum in JBJ| is over a discrete sequence of numbers 
tending to — oo. "Asymptotic" means that for any j 

r 

(7) f(x)-J2Mx) = o(\\x\\ j ) 

as 1 1 a; 1 1 — > oo. The number r occurring in (|SJ| is again called the order of the 
asymptotic series and the collection of functions satisfying © will be called poly- 
homogeneous symbols of order r and denoted by S r . 

For the sake of exposition in this introduction, we will continue to discuss the 
case where A is regular. 

We will show that if / is a polyhomogeneous symbol then the function p(N, f) 
given by Q is a polyhomogeneous symbol in N and its asymptotic expansion in 
powers of N is given by a formula similar to (JSJ with two key differences: 

(1) For symbols, an infinite number of differentiations occur on the right hand 
side of ©, i.e. the whole Todd operator must be applied. So JHJ must be 
understood as an asymptotic series, not as an equality. 

(2) The formula JSJ has to be corrected by adding a constant term C to the 
right hand side, a constant which is zero for the case of a polynomial. 

More precisely, we will prove: 

Theorem 1.1. Let A be a regular polytope whose vertices lie in Z n with in the 
interior of A. Let f S S r and N G Z + and let 

p(NJ):= Yl /(*)• 

Let p(t, h, f) be defined by 0) so that 

p(N,0,f)= f f(x)dx, 
Then p{N, f) — p(N, 0, /) is a symbol in N and has the asymptotic expansion 

p(N, f) - p{N, 0, /) ~ (toM (JA - ld\ p(N, h, f) 
where C is a constant. 



+ C 

h=0 



The constant C is of interest in its own right. It can be thought of as a "regu- 
larized" version of the difference 

(8) Yl /w - / ft*)**- 
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Of course there is no reason why either the sum or the integral in JHJ should 
converge. But we can "regularize" both as follows: Define the function (x) by 

{xf :=1 + ||*|| 2 . 

For s e C let 

f(x,a) :=f(x)(x} s . 

We will show that 

(9) C(s) := Y, /(^)- / ffas)dx, 

which is holomorphic for Re s << 0, has an analytic continuation to the entire 
complex plane and that the missing constant C on the right hand side of JHJ is 
exactly C(0). In particular, the constant C is independent of the particular polytope 
in question. 

This result is related to, and inspired by, a result of Friedlander-Guillemin |FG| 
on "Szego regularizations" of determinants of pseudodifferential operators. In their 
result, as in ours, there is a missing constant which also has to be computed by the 
above process of "zeta regularization" . 

Our result is somewhat insensitive to the mode of regularization. In fact, it can be 
generalized as follows: Define a "gauged symbol" |Ougau| to be a function f(x, s) £ 
C°°(IR™ x C) which depends holomorphically on s and for fixed s is a symbol of order 
Re s + r. For example, if / is a symbol, the function f(x){x) s introduced above is 
such a gauged symbol. We make a similar definition of "gauged polyhomogenous 
symbols". We will prove that if f(x, s) is a gauged polyhomogenous symbol with 

f(x) = f(x,0) 

then the function given by JHJ with this more general definition of f(x, s) again 
extends holomorphically from Re s << to the entire plane and C = C(0). 

The above results will be proved for the more general case of simple integral 
polytopes in §2. The proof is largely based on the Euler-MacLaurin formula with 
remainder as proved in KSW and motivated by an argument of Hardy on "Ra- 
manujan regularization" [Hardy]. Ramanujan's key idea was to use the classical 
Euler-MacLaurin formula in one variable to regularize (jHJ) by providing "counter 
terms" in passing to infinity in the difference between sum and integral in one di- 
mension. The origin of this method goes back to Euler's continuation of the zeta 
function past the pole at z = 1 and his introduction of what is known today as 
"Euler's constant". 

We would like to thank Yael Karshon and Michele Vergnc for helpful comments. 



2. An Ehrhart formula for simple polytopes. 

We continue with the notation of §1. So A C R" is a convex polytope whose 
vertices are in Z™ and such that the origin is in the interior of A. We assume in 
this section that A is simple, which means that for every vertex p there are exactly 
n edges emanating from p so they lie on rays 

(10) P + ta itP 

where the a^p, i = 1, . . . ,n form a basis of R™. 
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2.1. The Euler-MacLaurin formula for symbols. We want to apply the Euler- 
MacLaurin formula with remainder, |KSWj . In |KSWj one dealt with a weighted 
sum where points in the interior of the polytope are given weight 1, points on the 
relative interior of a facet are given weight w(x) := i, and, more generally, points 
in the relative interior of a face of codimension k are given weights w(x) := 
The weighted sum pi (N, /) is then defined as 

pi(NJ):= J2 

Theorem 3 of KSW gives an Euler-MacLaurin formula with remainder for weighted 
sums of symbols. 

More generally, |AW| consider the Euler-MacLaurin formula with remainder for 
more general weightings including the unweighted sum we considered in §1. We 
refer to equations (28) and (29) in |AW| for the definition of a general weighting, w, 
and we will denote the corresponding weighted sum here by p w (N, /). They stated 
their formula with remainder for smooth functions of compact support, but the 
passage from the case of smooth functions of compact support to that of symbols 
is exactly the same as in |KS\\ I . 

To formulate the Euler-MacLaurin formula with remainder we fix a vector £ £ K™ 
such that 

(11) cti iP ■ £ ^ V p and i. 

We then define 

if ot hp ■ £ > 
-a iyP if a hp • £ < 

a tt • £ 



(12) < p := 

(is) (-iy ■■= n 



=i l 



and 

(14) C P)t := j^ + X>4 p! *i>o|. 

There is a certain infinite order differential operator M (depending on the weight- 
ing) in the variables hi, ... , h rn with constant term 1 whose truncation of order k 
is given by the sum in equation (89) of KSW] (for weight |) and the sum in equa- 
tion (56) in |AW| (for general weights) such that for any symbol / of order r and 
k > n + r 

(15) Pw (N,f)-p(N,0J)= ((M [k] ) (J^j -ldjp(N,h,f) ^ +R k (f,N) 
where M' fe ' denotes the truncation of M at order k and 

/|-y|=nfc \ 

(16) R k (f,N) = ^(-lf / J2 ^,k D1 f dx 



| 7 |=fe 



where the <\f^ k are bounded piecewise smooth periodic functions. For an explicit 
expression for these functions sec KSW or AW . We will not need this here. If 
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the polytope is regular, and we use unweighted sums, the operator M is exactly 
the operator Todd of §1. 

Notice that M [k] -1 has no constant term, so ((M [fc] ) (^) - Id) p(N, h, f)\ h=0 
involves integration of derivatives of / over faces of the polytope. These faces are 
moving out to infinity as N — > oo. A derivative of a homogeneous summand in © 
is itself a homogenous function. So as N — > oo, the integral of this derivative over 
a face of the polytope is a homogenous function in N whose degree depends on the 
degree of this derivative and the dimension of the face. So the contributions of the 
polyhomogenous terms with sufficiently negative degree in 10 will be homogeneous 
terms of high negative degree in N in iJTSl). By the same token, each homogenous 
summand in (jfJJ will yield a finite number of terms to each order in (|15|l . So the 
left hand side of (|15|l will be polyhomogenous symbol in N. 

2.2. The key idea. We want to investigate the behavior of the remainder (|16ll as 
N -> oo. 

Since the origin is in the interior of the polytope, it is in the interior of the 
cone generated by the a, jP for any vertex p. This means that when we use the cti iP 
as a basis of R™ (based at the origin), the point p has strictly negative coefficients 
relative to this basis. When we flip the basis from a itP to af , the coefficients of 
those edges which are actually flipped become positive. 

Condition implies that the function x i— > £ • x has a unique minimum on the 
polytope, and that this minimum is achieved at a single vertex, q. At this vertex, 
no edges are flipped, and at any other vertex p ^ q at least one edge is flipped. So 
for p ^ q we have 



( 17 ) P = XXp a tp' 

i=i 

where at least one a^ p > 0. Then the cone C Pi t is contained in the half space 

\^X ^ ^ Xj Q^jp , Xi ^ G-i jP ^^ 

and so the p-th summand in i|16|) tends to zero as N — > oo for p ^ q. 

At the vertex q we have we have q — J^i a i,q a i,q = J2i a i,q a \ q w hh all the 
a.i :q < and so the cone C 9; t tends to the entire space R™ as t — > oo. Thus 
tends to 




It follows from l|15|l that this limiting value Ck is independent of k for k suffi- 
ciently large. Let us call this common limit C = C(f). It also follows from l|15|l 
that C(/) is independent of the choice of the polarizing vector £. 

We shall interpret this limiting value C using regularization in the next section, 
and we will find that C is also independent of the particular polytope we are 
expanding. 

If / is a polynomial, so that we choose k to be greater than the degree of /, we 
see from (|18|l that C = 0, as it must be from the classical Ehrhart theorem. 
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2.3. Application to polyhomogeneous symbols. Now suppose that / is a poly- 
homogenous symbol. We can apply the above to each summand in the asymptotic 
series ©• But notice that if we choose j sufficiently negative, the function 

r 

g(x) = g 3 (x) = f{x) ~J2M X ) 

i=j 

occurring on the left hand side of J7J will have the property that both 



Y 9(t) and J g(x)da 



are absolutely convergent. Furthermore, given any negative number m we can 
arrange, by choosing j sufficiently negative, that 



Pi 



and 



so 



p(N,0,g) 



g{x)dx = o{N m ) 



\p w (N,g)-p(N,0,g)}- 
So if we define 

r 

(19) C(/) = ^C(/,) 



E 

eel," 



g(x)dx 



o(N m ). 



gj{x)da 



for j sufficiently negative, then C(f) is independent of the choice of j. Furthermore, 
we see that if / is a polyhomogeneous symbol, we get an asymptotic expansion of 
the form 



p w (N,f)-p(N,0J)~ ( M 



d_ 

dh 



ld)p(N,h,fi 



+ C(f), 



h=0 



where each level in the asymptotic expansion in N involves only finitely many fi 
By abuse of language, we shall denote this equation as 



(20) p w (N,f)-p(N,0J) 



M 



d_ 

dh 



-Id \p(N,h,f) 



C(f). 



h=0 



2.4. Regularization. Suppose that we replace / by a gauged polyhomogenous 
symbol f(x,s) with f(x,0) € S r . Then the remainder term l|16fl applied to a 
summand in the asymptotic expansion of f s = f{-,s) is well defined if Re s < 
—r — n + k. Moreover, if p ^ q so that a^p > for some i the p-th summand on 
the right of (O is of order 0(jV Res + r +™-* ; ). 

At the unique vertex q where no edges are flipped g-th summand of 116|) differs 
from the integral 



(21) 



„ /h\=nk \ 

/ E 4>«, k Dlfi(x,s)\dx 
■ \ | 7 |=fc ' 



s 



V.W.GUILLEMIN, S. STERNBERG, AND J. WEITSMAN 



by a term of order 0(N Re s+i +n k ). Thus the gauged version of (121 III is 

d 



(22) Pw (NJ s )~p(N,0J s )= (M [fel -Id) (JL\p(N,hJ s ) 



dh / 

j r O(N Rcs+r+n ~ k ) 



+ C k (s) 



h=0 



for 



Re s < —n — r + k 

where f s = /(-, s) where Ck(s) is (|2()[) with / replaced by f s , and we have computed 
C(s) by going out to level k in the Euler-MacLaurin expansion. 

All the terms on the right of 1L'2I) are holomorphic on the half-plane Res < 
—71 — r + k. 

Letting k — > oo we conclude that on this half-plane we have, in the notation of 



d_ 

dh 



Id )p(N,h,f s 



C(s) 



fe=0 



(23) p w (NJ s )-p(N,0J s )~ 
where 

(24) C(s) - C fe (s) = o(N Rcs+r+n - k ). 

Since the Cfe(s) are holomorphic of the half-plane Re s < —n — r + k, it follows that 
C(s) is holomorphic on the whole plane. 

Moreover, in the asymptotic series on the right of all the terms are of order 
at most Re s + r + n. Hence for Re s < —r — n these terms tend to zero and we get 



C(s) = lim (p(N,f s )-p(N,0J s )) 

iV— >OG 

= / f(x,s)dx, 



and both the sum and the integral converge absolutely. So if we set s = we obtain 

d 



(25) p w (N, f) - p(N, 0, /) ~ ( M ( ^ ) - Id ) p(N, h, f) 



C 



h=0 



where f(x) — f(x,0) and C — C(0). So we can think of C as a "regularization" of 
To summarize: We have proved 

Theorem 2.1. Let A be a simple polytope whose vertices lie in Z™ with in the 
interior of A. Let f S S r and N S Z+ and let 

p w (N,f):= J2 W (W)- 

Let p(t, h, f) be defined by ^ so that 

p(N,0,f)= f f(x)dx, 

JN-A 

Thenp(N,f) — p(N, 0,/) is a symbol in N and has the asymptotic expansion 

d 



p w (N,f)-p(N,0J) 



M 



Oh 



Id )p(N,h,f) 



c 



h=0 
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where C is a constant. Furthermore, if f{x, s) is a gauged polyhomogenous symbol 
with f(x,0) = f{x) then C = C(0) where C(s) is the entire function given by I2tyl 
and \2J$ - For Re s < —r — n 

C ( s ) = E s ) - / f(x,s)dx. 

Hence C(s) and in particular C — C(0) is independent of the polytope. 

Remarks. 1. In the course of the discussion we have proved a similar theorem 
with polyhomogeneous symbols replaced by symbols and gauged polyhomogenous 
symbols replaced by gauged symbols. 

2. Since the initial posting of this paper we have received the interesting paper 
[MP] . 
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